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This  study  experimentally  and  numerically  investigates  the  effects  of  the  four  different  types  of  serpentine 
flow-field  geometry  on  the  performance  of  a  direct  methanol  fuel  cells  (DMFCs).  To  elucidate  the  effect 
of  different  channel/rib  aspect  ratios,  the  through/in-plane  transport  of  methanol  and  air,  the  current 
density  distribution,  the  anode  and  cathode  polarization,  the  impedance  and  the  cathode  pressure  drop 
are  numerically  and/or  experimentally  observed.  The  sub-rib  convection  significantly  affects  the  cell 
performance  with  the  serpentine  flow  fields.  Thus,  a  narrow  rib  width  and  a  suitable  channel/rib  aspect 
ratio  are  needed  to  increase  the  cell  performance  of  serpentine  DMFCs.  This  flow  channel,  for  example, 
an  F2-type  flow  channel  (channel  width:  1 .0  mm,  rib  width:  0.5  mm),  leads  to  formation  of  a  narrower  rib 
chain  as  well  as  a  higher  current  density  distribution  in  the  rib,  and  thereby  it  brings  about  improvement 
of  the  two-phase  flow  in  DMFCs,  which  is  supported  through  the  high  performance  difference  in  the  high 
current  density  regions  and  the  long-term  stability  test  of  700  h. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  methanol  fuel  cells  (DMFCs)  are  attractive  as  a  promising 
power  source  for  portable  multifunctional  devices  and  light-duty 
vehicles  because  of  easy  storage  of  the  fuel,  light  weight,  small 
system  size  and  high  energy  density  [  1  ].  The  performance  and  dura¬ 
bility  of  DMFCs  depend  on  many  factors,  including  the  operating 
conditions,  mass  transport  phenomena  inside  the  cell,  kinetics  of 
the  electrochemical  reaction,  membrane  electrode  assembly  (MEA) 
and  flow  channel  geometry  [2]. 

The  flow  field  channel,  which  is  integrated  into  bipolar  plates, 
is  also  an  important  component  because  it  serves  as  a  fuel  or  air 
distributor  as  well  as  a  current  collector.  In  general,  the  flow  field 
channels  influence  the  fuel  and  air  distributions,  the  concentra¬ 
tion  polarization  loss,  the  electrical  current  conduction,  the  heat 
conduction,  etc.  [3-9].  The  flow  field  channel  is  especially  impor¬ 
tant  for  water  management  inside  the  fuel  cells  since  the  products 
(CO2  and  H20)  are  removed  through  flow  field  channels.  There¬ 
fore,  an  appropriate  flow  channel  design  is  necessary  for  achieving 
high  performance  and  durability  of  the  cell  as  it  resolves  mass 
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transport-related  problems  and  improves  the  water  management 
in  the  channel,  the  gas  diffusion  layer  (GDL)  and  the  catalyst  layer 
(CL)  [10]. 

Many  researchers  have  studied  the  effects  of  flow  channel 
geometry  on  the  homogeneous  distribution  of  the  reactants  and 
on  cell  performance.  Shimpalee  et  al.  [11]  investigated  the  effects 
of  the  channel  and  rib  cross-sectional  areas  on  the  reactant  distri¬ 
bution  and  cell  performance  with  serpentine  flow  field  patterns  and 
showed  that  narrower  channels  with  wider  rib  spacing  resulted  in 
higher  performance.  Yoon  et  al.  [12]  reported  that  a  narrower  rib 
width  led  to  improved  performance  of  proton  exchange  membrane 
fuel  cells  (PEMFCs)  because  the  diffusion  of  the  gas  appeared  to 
be  a  more  important  factor  than  the  electrical  conduction  for  bet¬ 
ter  cell  performance.  Wang  et  al.  [13,14]  investigated  the  effects 
of  the  number  of  flow  channel  bends,  the  number  of  serpentine 
loops,  and  the  flow  channel  width  ratio  on  the  cell  performance  of 
PEMFCs  with  serpentine  flow  fields.  They  particularly  probed  the 
effect  of  the  sub-rib  convection  on  cell  performances  under  var¬ 
ious  channel  aspect  ratios.  According  to  their  results,  the  sub-rib 
convection  of  the  serpentine  flow  field  significantly  affects  the  reac¬ 
tant  transport  efficiency  through  the  GDLs  to  CLs,  and  the  removal 
rate  of  liquid  water  accumulated  underneath  the  ribs.  Manso  et  al. 
[15]  numerically  investigated  the  influence  of  the  channel  cross- 
section  aspect  ratio  (defined  as  the  ratio  height/width)  on  the 
performance  of  a  PEMFC  with  a  serpentine  flow  field  design,  and 
analyzed  the  local  current  density  distributions  in  the  membrane, 
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Nomenclature 

a  water  activity  or  effective  catalyst  area  per  unit  of 

the  total  volume  (m2  m-3) 

A  area(m2) 

C  molar  concentration  (mol  m-3 ) 

Dj  mass  diffusivity  of  species  i  (m2  s-1 ) 

EW  equivalent  weight  of  a  dry  membrane  (kg  mol-1 ) 

F  Faraday  constant  (96487  C  mol-1 ) 

10  exchange  current  density  (A  m-2 ) 

I  operating  current  density  (A  m-2 ) 

I xover  crossover  current  density  (A  nrr 2 ) 
j  transfer  current  density  (A  m-3 ) 

jl  diffusive  mass  flux  of  phase  i  (kg(m2  s)-1 ) 

J  Leverett  function 

kr  relative  permeability 

I<  hydraulic  permeability  (m2) 

M  molecular  weight  (kg  mol-1) 

m  mass  fraction 

n  number  of  electrons  in  the  electrochemical  reaction 

or  the  diffusivity  correction  factor 
nc  catalyst  coverage  coefficient 

nd  electroosmotic  drag  coefficient 

n  xover  crossover  molar  flux  (mol(m2  s)-1 ) 

P  pressure  (Pa) 

Pc  capillary  pressure  (Pa) 

Ru  universal  gas  constant  (8.314 J  (mol  I<)-1 ) 

s  stoichiometry  coefficient  in  the  electrochemical 

reaction  or  liquid  saturation 
S  source  term  in  the  transport  equation 

t  time  (s) 

T  temperature  (K) 

u  fluid  velocity  and  superficial  velocity  in  a  porous 

medium  (ms-1) 

U0  thermodynamic  equilibrium  potential  (V) 

V  volume  (m3) 

Vceu  cell  potential  (V) 

Greek  letters 

a  transfer  coefficient 

8i  thickness  of  component  i 

e  volume  fraction  of  the  gaseous  phase  in  the  porous 

region 

8mc  volume  fraction  of  the  ionomer  phase  in  the  CL 

y  advection  correction  factor 

X  membrane  water  content  (molH2o(molso3-  )_1 ) 

Xa  relative  mobility  of  phase  a 

0  phase  potential  (V) 

r 7  overpotential  (V) 

0  contact  angle  (°) 

11  viscosity  (kg(m  s)-1 ) 

p  density  (kg  m-3) 

pmem  dry  membrane  density  (kg  m-3 ) 

v  kinematic  viscosity  (m2  s-1 ) 

a  surface  tension  (Nm-1)  or  electronic  conductivity 

(Sm-1) 

r  viscous  shear  stress,  N  m-2,  tortuosity 

k  ionic  conductivity  (S  m-1 ) 

§  stoichiometry  flow  ratio 

Superscripts 
e  electrolyte 

eff  effective  value  in  the  porous  region 

mem  membrane 


8  gas 

l  liquid 

ref  reference  value 

s  solid 

sat  saturation  value 

xover  crossover 

Subscripts 
a  anode 

avg  average  value 

BP  bipolar  plate 

c  cathode 

CL  catalyst  layer 

e  electrolyte 

GC  gas  channel 

GDL  gas  diffusion  layer 

MeOH  methanol 

i  species  index 

in  channel  inlet 

m  mass  equation 

mem  membrane 

02  oxygen 

ref  reference  value 

t  total 

s  solid 

sat  saturation  value 

u  momentum  equation 

w  water 

0  potential  equation 

0  standard  condition,  298.1 5  K  and  1 01 .3  kPa  ( 1  atm) 


velocity  distributions,  liquid  water  saturation  and  hydrogen  and 
oxygen  concentrations  to  understand  the  channel  cross-section 
aspect  ratio  effect.  These  studies  demonstrated  that  flow  channel 
geometry  influences  the  performance  of  PEMFCs  in  terms  of  fuel 
distribution,  mass  transport  and  water  removal. 

In  DMFCs,  the  flow  channel  is  a  critical  component  for  mini¬ 
mizing  the  mass  transport-related  issue  raised  by  the  two-phase 
flow  of  liquid  (CH3OH  and  H20)  and  gas  (C02  and  air)  and  for 
removing  a  large  amount  of  liquid  water  on  the  cathode  side,  which 
results  from  product  water  by  oxygen  reduction  reaction  (ORR)  and 
electro-osmotically  dragged  water  from  the  anode.  Yang  et  al.  [16] 
experimentally  examined  the  effect  of  the  anode  flow  field  pattern 
on  the  performance  of  a  DMFC  and  reported  that  both  the  open 
ratio  and  the  flow  channel  length  had  a  significant  influence  on  the 
cell  performance  and  the  pressure  drop.  Despite  these  important 
facts,  there  have  been  few  papers  that  dealt  with  the  effect  of  the 
flow  channel  geometry  on  the  performance  of  DMFCs. 

Thus,  this  paper  numerically  and  experimentally  investigates 
the  effects  of  flow  channel  geometry  on  the  performance  of  DMFCs 
using  a  single-path  serpentine  pattern  with  different  channel  and 
rib  widths  and  channel  lengths.  The  through/in-plane  transport  of 
methanol  and  air  in  DMFCs,  the  current  density  distribution  and 
cathode  pressure  drop  are  numerically  analyzed  to  examine  the 
effect  of  the  flow  channel  cross-sectional  area  on  the  cell  per¬ 
formance  for  single  serpentine  DMFCs.  In  addition,  a  single  cell 
test,  an  anode/cathode  polarization  test,  an  impedance  analysis  and 
measurement  of  the  cathode  pressure  drop  are  experimentally  con¬ 
ducted  to  evaluate  the  characteristics  of  the  different  flow  channel 
geometries  in  DMFCs. 
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2.  Numerical  model 

2.1.  Multi-dimensional,  multi-phase  fuel  cell  model  and 
assumptions 

The  two-phase  DMFC  model  used  in  this  study  is  based  on  the 
multi-phase  mixture  (M2)  formulation  of  Wang  and  Cheng  [17]. 
The  present  numerical  model  is  developed  from  the  PEFC  model 
developed  earlier  [18-22]  and  modified  for  use  in  a  DMFC  system. 
Utilizing  the  M2  formulation  to  describe  two-phase  transport  in 
DMFCs,  the  present  model  makes  the  following  assumptions: 

(1)  Laminar  flow  due  to  small  flow  velocities. 

(2)  Ideal  gas  mixtures  in  the  gas-phase  region. 

(3)  Isotropic  porous  media  (GDLs  and  CLs). 

(4)  Negligible  effect  of  immobile  liquid  saturation  in  the  porous 
diffusion  medium. 

The  last  assumption  is  based  on  the  study  of  Ju  [19]  in  which 
the  effect  of  immobile  liquid  saturation  is  diminished  when  the 
diffusion  media  faces  a  high  liquid  saturation  level. 

2.2.  Conservation  equations 


In  Eq.  (8),  p£  is  the  gas  mixture  density  described  by  the  ideal 
gas  law;  it  varies  with  the  species  composition  (the  mass  fractions, 
{mf  }).  That  is, 


^  (rut)  Y.mf/M'- 

(9) 

Mixture  velocity: 

-  plul  +  pgug 

u  =  - - - - . 

p 

(10) 

Species  mass  fraction: 

plsm\  +  p£(l  -  s)mf 

m,  = - - - - . 

P 

(11) 

Relative  permeability: 

k[  =  s4, 

(12a) 

-?6 

ii 

i 

> 

(12b) 

Kinematic  viscosity: 

7 

+ 

ii 

(13) 

Under  the  above  assumptions,  the  two-phase  DMFC  model 
is  governed  by  the  conservation  of  mass,  momentum,  chemical 
species  and  charge. 

Mass  conservation: 

V  •  (pu)  =  Sm  (1) 

Momentum  conservation: 

Flow  channels  (Navier-Stokes  equations). 

V  •  (puu)  =  —  Vp  +  V  •  r  (2) 

Porous  media  (Darcy’s  equations). 

pu  =  "Vp  (3) 

Species  conservation: 

Flow  channels  and  porous  media. 

V  ■  ( YipniiU )  =  V  ■  [yOgDf'e-flrV(mf )]  +  V  •  [plDl:effV(m\)] 

+  V-[(mf  -m\tf]+St  (4) 

Water  transport  in  the  membrane. 

(nmem  \  /  1  \  /  vunem  \ 

WD"em  W  -  V  •  ("4 F )  Mw  +  V  •  (~VP'J  =  0 

(5) 

Charge  conservation: 

Proton  transport :  V  •  (/ce^V0e)  +  S&  =  0.  (6) 

Electron  transport :  V  •  (ore^V<2>s)  -  S&  =  0  (7) 

2.3.  Two-phase  transport  properties 

Generally,  two-phase  mixture  properties  are  expressed  as  func¬ 
tions  of  s,  the  volume  fraction  of  pores  occupied  by  the  liquid  phase, 
and  (1  -s),  the  volume  fraction  of  pores  occupied  by  the  gaseous 
phase.  The  different  mixture  properties  are  defined  as  follows  [17]. 
Mixture  density: 


(8) 


In  Eq.  (13),  v§  is  the  kinematic  viscosity  of  the  gas  mixture  that 
varies  with  the  gas  composition  [23]: 


ps 


Xi/ii 


,x,<V 

i=i  j  y 


where 


0ij  = 


i 

V8 


and 


jii  (N  s  m  2)  = 


fiw  =  0.00584  xlO"6!"1  29 

AIn2  =  0.237  x  10-6T0-76  ,  T  in  K. 

fio2  =  0.246  x  10-6T°-78 


(14) 


Relative  mobility: 

kl  =  —j  v,  (15a) 

vl 

X8  =  \-kl.  (15b) 


The  species  diffusivity  in  the  gas  mixture,  Df,  in  the  first  term 
in  the  right-hand  side  of  Eq.  (4)  is  defined  as  follows  in  order  that 
the  cumulative  interspecies  diffusion  within  the  gaseous  phase  is 
equal  to  zero  [23]. 


D?  = 


1  -X[ 


£7  VDy 


j±i 


where  Djj  = 


1.013  x  10“7  -T1-75 

1/3^  1/3  r 


p-ur+xjn 

Xw  =  12.7,  xn2  =  17-9,  Xo2  =  16.6 


1  1 

Wi  +  W] 


1/2 


(16) 


Note  that  the  gaseous  diffusive  transport  for  a  porous  medium 
can  be  controlled  by  the  Knudsen  diffusion  effect  due  to  molecule- 
to-wall  collision  as  well  as  the  molecular  diffusion  which  is  caused 
by  inter-molecular  collision,  as  described  in  Eq.  (16).  The  Knud¬ 
sen  diffusion  coefficient  can  be  computed  according  to  the  kinetic 
theory  of  gases  as  follows. 


K_  2  fSRuT\ 

1  3  V  7TM;  J 


1/2 


p  =  pls  +  pg(  1  -s). 


(17) 
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The  effective  diffusivity  of  species  i  in  the  gas  mixture  is  then 
obtained  by  combining  both  the  molecular  and  Knudsen  diffusion 
effects  with  the  effects  of  the  porosity  and  tortuosity  of  the  porous 
medium  using  Bruggeman’s  correlation. 

of*  =  Ws)r(±  +  ^)  ns) 


The  effective  diffusivity  of  species  i  in  the  liquid  phase  is 
obtained  as, 

Dfff  =  [es]rD|  (19) 

Assuming  the  vapor-liquid  equilibrium  of  methanol  on  the 
anode  side  and  invoking  Henry’s  law,  the  transportation  equation 
for  methanol  can  be  written  as, 


V  ■  (YMeOH pmMeOHU)  — 

Pll<H  ) 

VmL 

j  jxover 

SMeOH  =  -  g p  -  —Qp~  Anode  catalyst  layer 


kH 


QS,eff 

MeOH  \  v(  / 

^  I  Vk"lMeOHJ 


+  SMeOH 


(20) 


where  Henry’s  constant,  /<#,  is  defined  as  the  ratio  between 
the  methanol  concentration  in  the  liquid  and  gas  phases  and  is 
expressed  as  a  function  of  temperature,  i.e., 


ri 

]<H  =  ^ MeOH  =  kOe-[(AHsol/R)((l/T)-(l/T0))]RT ,  (21) 

CjvieOH 

where  k ^  refers  to  Henry’s  constant  at  the  reference  state  of 
T0  =  298.1 5  K  and  A Hso/  refers  to  the  enthalpy  of  the  solution.  It 
is  seen  that  the  effective  liquid  and  gas  diffusivities  in  the  two- 
phase  region  are  functions  of  both  of  the  porosity,  £,  and  the  liquid 
saturation,  s. 

Assuming  complete  consumption  of  methanol  at  the  cathode 
catalyst  layer  after  it  crosses  over  the  membrane  from  the  anode  to 
cathode,  the  methanol  flux  driven  by  electro-osmotic  drag  (EOD) 
and  diffusion  can  be  written  as: 


J(s)  = 


I1' 

u 


417(1  -s)- 2.120(1  -s)2  + 1.263(1  - 
417s  -  2.120s2  +  1.263s3 


sf 


if  0C  <  90° 
if  0C  >  90° 


(26) 


The  term  in  the  left-hand  side  of  the  species  equation  represents 
the  advective  term,  in  which  the  advection  correction  factor,  yit  is 
given  by  Wang  and  Cheng  [17]. 


Yi 


s  = 


p(Xlml  +  X8m8) 

(. splm\  +  (1  -  s)psm :?) 

The  liquid  saturation,  s,  can  be  calculated  via, 

r  _  rS 
'-w  ^ w,sat 


Cl  -  C8 

'-w  ^ w,sat 


(27) 


(28) 


where  liquid  water  concentration,  Clw,  can  be  calculated  by  assum¬ 
ing  the  vapor-liquid  equilibrium  of  methanol  as  follows: 


rl  _  QweOH 

^oh-  [s  +  ((l-s)//<H)] 

/  _  ^MeOH^MeOH 

Me0H~  [(S  +  ((l-S)/kH))p'] 
,  _  (1  -  P^MeOH^P1 

-  M 


The  transport  properties  of  electrolytes  are  correlated  with  the 
water  content  of  the  membrane,  X,  which,  in  turn,  is  a  function  of 
the  water  activity,  a,  as  follows  [24]. 


cJaj 

Psat 


(30) 


/  =  0.043  +  17.81a-  39.85a2  +  36.0a3  for  0<a<l,_1. 

\\'  =  22  (31} 


The  electro-osmotic  drag  coefficient,  nd ,  the  proton  conductiv¬ 
ity  in  the  membrane,  /c,  and  the  water  diffusion  coefficient  in  the 
membrane,  D™em,  are  given  by  Springer  et  al.  [24]. 

2-5A 

^d,w  —  22  (32) 


k  =  (0.5139A  -  0.326)exp 


1268 


(33) 


r>  mem 
uw 


'  2.692661843  x  10“10 

{0.87(3  -2.)+  2.95(7.  -  2))  X  10“10  •  e(7.9728-(2416/T)) 

{2.95(4  -  7.)  +  1 .642454(7.  -  3))  x  10“10  •  e(7-9728-(24ie/T)) 

.  (2.563  -0.337. +  0.0264A2  -0.00067U3)  x  10“10  •  e(7-9728-(2416/7')) 


for  X  <2 
for  2  <  X  <  3 
for  3  <  A  <  4 
for  4  >  X 


(34) 


6  F 


=  V, 


(^nd,MeOH^  + 


nmem 

uMeOH 


Cl  I 

MeOH  |  cata 


(22) 


where  nd  Me0H  denotes  the  EOD  coefficient  for  methanol  which  is  a 
function  of  EOD  coefficient  of  water,  and  the  ratio  of  the  methanol 
to  water  concentrations  in  the  anode  CL.  That  is, 


nd,MeOH  =  nd,w 


Cl 

^ MeOH 


aCL 


(23) 


The  diffusive  mass  flux  of  the  liquid  phase,/,  which  is  shown 
in  the  third  term  in  the  right-hand  side  of  Eq.  (4),  is  expressed  as  a 
function  of  the  capillary  pressure,  Pc,  such  that  [17]: 


j1  =  plul  -  x'pu  =  ^X'XgVPc 


(24) 


In  Eq.  (24),  the  correlations  of  the  capillary  pressure,  Pc,  and  the 
Leverett  function  J(s),  are  found  as  below. 


Pc 


ps  -Pl  =  cr cose 


1/2 

J(S), 


(25) 


2.4.  Electrochemical  reactions  and  source  terms 

In  the  above  conservation  equations,  namely,  (l)-(7),  Sm,  Sz 
and  S&,  which  are  listed  in  Table  1,  denote  the  corresponding 
source/sink  terms  for  mass,  species  and  charge.  These  source  terms 
result  from  the  electrochemical  reactions:  the  methanol  oxidation 
reaction  ( MOR)  in  the  anode  CL  and  the  ORR  in  the  cathode  CL.  These 
electrochemical  reactions  are  represented  by  kinetic  expressions  as 
follows. 

MOR  in  the  anode  CL  [25-27]: 

J  =  CMeOH  |cataexP((«qF/R7-)^) 

'-MeOH  |  cata  +  $  exP ((aaF/RT)r]a) 

ORR  in  the  cathode  CL: 

/  r  \  3/4 

k  =  -(1  -s)ncaifc  exp  (-^yF^)  (44) 

where 

jc=j+jxover  (45) 


Table  1 

Two-phase  steady-state  DMFC  model:  source/sink  terms. 


36 


Y.-C.  Park  et  al.  /  Journal  of  Power  Sources  205  (2012)  32-47 


r^ 

ro 


+ 


> 


+  m 

O 

II 

O  + 
Xx 

+ 

X  + 
O  JN 


Q 

"qj 

> 


W 


c  2 

rS  ^ 

■S  3 
S  £ 
g* 

a 

^  O 


rc  -S5 

S  tJ 

X  3 
O  T3 
—  QJ 
O  , 
C  C 
03  CU 
bJO 

*  2  X 

S  O 


In  Eq.  (44),  (1  -s)nc  approximates  the  effect  of  the  decreasing 
number  of  electrochemically  active  catalyst  sites  due  to  the  pres¬ 
ence  of  water  in  the  CLs. 

The  surface  overpotential  in  Eqs.  (43)  and  (44)  is  defined  by: 
MOR  in  the  anode  CL: 

rj  =  (ps-(pe  (46) 

ORR  in  the  cathode  CL: 

r)  =  (ps-(Pe-  Uo  (47) 

In  Eq.  (47),  the  thermodynamic  equilibrium  potential  between 
the  cathode  and  anode  can  be  written  as  follows. 

Uo  =  1.23  -  0.9  x  10-3(T  -  298.15)  (48) 

The  temperature  dependence  of  the  MOR  kinetics  [28]  and  ORR 


kinetics  [29]  are  expressed  as  follows. 

di'fjT)  =  aifa(353K)  •  exp 

r  35570  / 

1 

1  V 

(49a) 

L  R  V 

T 

353.15 /_ 

air0ef  T)  =  00353/0  •  exp 

’  73200  ( 

1 

1  T 

(49b) 

\  R  V 

T 

353.15  /  J 

2.5.  Boundary  conditions 

The  inlet  velocities  in  the  anode  and  cathode  gas  channels  can 
be  expressed  as  functions  of  the  anode  and  cathode  stoichiometric 
ratios  ($a  and  §c),  respectively,  the  DMFC  operating  current  density 
(/),  the  cross-sectional  areas  of  the  anode  and  cathode  gas  channels 
(. Aa  and  Ac),  respectively,  and  the  inlet  concentration  of  methanol 
and  oxygen,  which  are  functions  of  the  anode/cathode  inlet  pres¬ 
sure,  temperature,  and  humidity. 


$a(I/6F) Amem 

^MeOH^a 


and 


Uin,c  ~ 


&(//  4  F)Amem 
Co2Ac 


(50) 


In  addition,  the  isothermal  boundary  condition  (60 °C  in  this 
study)  and  the  constant  current  boundary  condition  are  applied 
to  the  bipolar  plate’s  outer  surfaces,  whereas  no-slip  and  imper¬ 
meable  velocity  condition  and  no-flux  condition  are  applied  to  the 
outer  faces  for  the  velocity  and  species  calculations. 


2.6.  Flow  channel  geometry  and  numerical  procedures 


The  DMFC  model  described  above  was  implemented  with  a 
commercial  computational  fluid  dynamics  (CFD)  package,  STAR- 
CD,  based  on  its  user-coding  capability  [30]. 

To  investigate  the  effects  of  the  flow  channel  geometry  on  the 
performance  of  DMFCs,  single-path  serpentine  patterns  with  an 
active  area  of  9  cm2  (3  cm  x  3  cm)  was  conducted  for  the  numerical 
analyses.  The  flow  channel  geometries  were  designed  with  differ¬ 
ent  channel  and  rib  widths,  channel  lengths,  cross-sectional  area 
ratios  and  number  of  U-turns,  as  shown  in  Fig.  1 .  Numerical  meshes 
were  generated  for  four  different  geometries  with  a  single-path  ser¬ 
pentine  flow  field.  The  flow  channels  had  channel  widths  of  0.7, 1.0 
and  1.5  mm,  rib  widths  of  0.5  and  1.0  mm,  and  channel  lengths  of 
450, 570  and  750  mm;  all  of  the  channels  had  a  depth  of  1 .0  mm.  The 
detailed  physical  geometries  of  the  single-path  serpentine  are  sum¬ 
marized  in  Table  2.  The  flow  field  designs  with  channel/rib  widths 
of0.7/0.5mm,  1. 0/0.5  mm,  1. 5/0.5  mm  and  1. 0/1.0  mm  are  denoted 
as  F1-,  F2-,  F3-  and  F4-type  channel,  respectively. 

Fig.  2  shows  the  mesh  configuration  of  the  F2-type  geometry.  In 
the  present  work,  the  computational  domain  of  a  DMFC  consists  of 
seven  sub  regions,  namely,  flow  channels,  backing  layers  (BL)  and 
CLs  in  both  anode  and  cathode  sides,  and  a  membrane.  Roughly 
1.2  million  cells  were  applied  for  each  of  the  four  geometries  and 
it  required  16  s  per  iteration  on  a  8-PC  cluster  (Intel®  Core™  i7 
processor  2.53  GHz). 
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Table  2 

Channel  name  and  physical  geometries  of  the  flow  field  patterns  with  different  channel  and  rib  widths. 
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Channel  name 

Channel/rib  width  (mm) 

Channel  length  (mm) 

Channel  area  (mm2) 

Rib  area  (mm2) 

Percentages  of  rib 
area  in  active 
electrode  area  (%) 

FI 

0.7/0.5 

750 

525 

360 

40.7 

F2 

1. 0/0.5 

570 

570 

270 

32.1 

F3 

1. 5/0.5 

450 

675 

210 

23.7 

F4 

1. 0/1.0 

450 

450 

420 

48.3 

Fig.  1.  Flow  channel  geometries:  (a)  F1-,  (b)  F2-,  (c)  F3-  and  (d)  F4-type. 


Fig.  2.  Mesh  configuration  of  the  F2-type  geometry. 
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3.  Experimental 

3.1.  MEA  preparation 

The  MEA  with  an  active  electrode  area  of  9.0  cm2  was  fab¬ 
ricated  using  commercial  Nafion  115  as  the  membrane,  PtRu/C 
(HISPEC  12100,  Johnson  Matthey  Fuel  Cells,  UK)  as  the  anode  cat¬ 
alyst,  and  Pt/C  (HISPEC  13100,  Johnson  Matthey  Fuel  Cells,  UK)  as 
the  cathode  catalyst  [31].  The  catalysts,  Nafion  solution  (10wt%, 
DuPont,  USA),  distilled  water,  isopropyl  alcohol,  and  1 -propanol 
were  mixed  to  make  the  electrode  catalyst  slurry,  and  the  mix¬ 
ture  was  homogenized  with  an  ultrasonic  processor  (Hielscher, 
UP100H,  Germany).  The  mixture  was  then  stirred  with  a  magnetic 
stirring  system  to  prepare  the  appropriate  catalyst  slurry.  The  cata¬ 
lyst  slurry  was  painted  onto  the  GDL  using  a  brushing  method.  The 
GDL  for  the  anode  and  the  cathode  were  carbon  paper  (TGP-H-060, 
Toray,  Japan)  that  was  treated  with  5  wt%  PTFE  and  a  carbon  sheet 
(GDL  25BC,  SGL  Carbon,  Germany),  respectively.  The  loading  of  the 
catalysts  was  2  mg  cm-2  PtRu  for  the  anode  and  2  mg  cm-2  Pt  for 
the  cathode.  The  electrode  layers  and  Nafion  membrane  were  hot 
pressed  using  a  laboratory  press  (Carver,  Model  M,  USA).  The  hot 
pressing  of  the  MEAs  was  performed  at  a  temperature  of  150  and  a 
pressing  pressure  of  50  kg  cm-2  for  60  s. 


3.2.  Characteristics  analysis  of  flow  channel  geometry 


To  investigate  the  effects  of  the  flow  channel  geometry  on  the 
performance  of  DMFCs,  the  single  cells  with  the  same  active  elec¬ 
trode  area  of  9.0  cm2,  GDL,  CL,  and  proton  exchange  membrane  was 
tested  using  an  electrochemical  test  system  (Won-A  Tech.,  Korea). 
A  solution  of  1  M  methanol  solution  was  fed  into  the  anode  at  a 
flow  rate  of  A  =  2.5  and  air  was  supplied  to  the  cathode  at  a  flow 
rate  of  A  =  3  under  ambient  pressure.  The  anode  polarization  of  the 
cells  was  also  measured  under  the  same  conditions  as  the  single 
cell  test;  however,  to  create  a  dynamic  hydrogen  electrode  (DHE) 
in  this  case,  hydrogen  (H2 )  gas  was  supplied  to  the  cathode  at  a  flow 
rate  of  10  ml  min-1. 

The  ohmic  resistance  and  electrochemical  impedance  spec¬ 
tra  of  the  single  cells  during  the  electrochemical  reaction  can 
be  observed  using  an  electrochemical  analysis  instrument  (Zah- 
ner,  IM6&IM6eX,  Germany).  Nyquist  plots  of  the  impedance  were 
recorded  under  galvanostatic  conditions  of  100  and  400  mAcm-2. 
All  anode  impedances  were  measured  between  the  anode  and  a 
DHE  in  the  complete  fuel  cell,  whereas  the  impedance  of  a  complete 
cell  was  obtained  using  a  DMFC  single  cell  system  (methanol/air). 
The  frequency  range  varied  from  3  mHz  to  1 0  kHz  and  the  ampli¬ 
tude  of  the  AC  voltage  was  1 0  mV. 

The  analytical  pressure  drop  of  the  different  flow  channel 
geometries  was  calculated  by  the  following  Darcy-Weisbach  equa¬ 
tion  [32], 


(51) 


Here,  A P  is  the  pressure  loss  due  to  friction, /is  the  friction  factor, 
L  is  the  channel  length  (m),  Dh  is  the  hydraulic  diameter  (m),  p  is 
the  fluid  density  (air  =  1 .23  kg  m-3 ),  and  v  is  the  velocity  of  the  flow 
(m  s-1 ).  For  the  laminar  flow  of  Re  <  2000,  the  friction  factor  (J)  is 
defined  as  Eq.  (52), 

'-is  (52> 

where  Re  is  the  Reynolds  number. 

Re  is  defined  as  the  following  Eq.  (53), 


Re  = 


PvDh 

li 


(53) 


where  /z  is  the  viscosity  of  the  fluid  (kgm-1s-1) 
(air  =1.81  x  10-5  kgm-1  s-1). 

For  a  rectangular  channel  with  wc  as  the  width  and  dc  as  the 
depth,  hydraulic  diameter  (Dh)  is  defined  as  Eq.  (54): 


Dh  = 


2  wc.dc 
wc  +  dc 


(54) 


Dh  and  v  as  well  as  Re  values  for  the  calculation  of  the  analyt¬ 
ical  pressure  drop  of  the  different  flow  channel  geometries  are 
represented  in  Table  3.  The  experimental  pressure  drop  in  the 
flow  channels  was  measured  at  each  current  density  from  100  to 
500  mA  cm-2  under  a  flow  rate  of  A  =  3  using  a  nanometer  that  was 
installed  at  the  inlet  of  the  cathode. 


4.  Results  and  discussion 

In  this  study,  the  presented  two-phase  DMFC  model  was  applied 
to  a  single-path  serpentine  channel  cell,  as  shown  in  Fig.  1  with  its 
mesh  configuration.  The  methanol  solution  and  air  were  fed  into 
the  anode  and  cathode  inlets,  respectively,  under  the  counter-flow 
configuration.  The  individual  cell  properties  and  operating  condi¬ 
tions  are  listed  in  Table  4,  and  the  relevant  physical  and  transport 
properties  are  given  in  Table  5.  Using  the  cell  configuration,  the  cell 
properties,  the  operating  conditions  and  the  physical  properties, 
a  numerical  study  was  conducted  to  investigate  the  effects  of  the 
flow  channel  cross-sectional  area  on  the  performance  of  serpentine 
DMFCs. 

In  general,  DMFCs  supply  methanol  and  air  through  the  flow 
field  channel  to  the  GDL  and  the  CL,  and  the  produced  C02  gas 
or  H20  from  the  cell  is  reversely  transported  out.  This  two-phase 
flow  behavior  of  DMFCs  significantly  affects  the  mass  transport  of 
reactants  (methanol  and  air)  and  the  removal  of  products  (C02  gas 
or  H20)  in  the  anode  and  cathode  side,  and  thereby  it  determines 
the  cell  performance  of  DMFCs  [36].  During  the  fuel  cell  operation, 
the  methanol  and  air  transport  processes  in  DMFCs  was  confirmed 
through  the  methanol  concentration  in  the  anode  and  the  air  con¬ 
centration  in  the  cathode. 

4.1.  Effect  of  the  flow  channel  geometry  on  performance  of  DMFCs 

Fig.  3  shows  the  experimentally  measured  and  numerically  sim¬ 
ulated  polarization  curves  of  the  single  cells  using  the  different 
flow  channel  geometries  listed  in  Table  2.  Both  the  numerical  and 
experimental  results  revealed  that  these  cells  exhibit  the  differ¬ 
ent  performances  with  the  flow  channel  geometries  employed  in 
each  cell.  According  to  the  experimental  results,  F2-type  flow  chan¬ 
nel  is  superior  to  other  type  flow  channels.  Based  on  the  same  rib 
width  of  0.5  mm,  the  F2-type  flow  channel  (channel  width:  1 .0  mm, 
rib  width:  0.5  mm)  had  the  highest  maximum  current  density  of 
670 mAcm-2  at  0.2V  and  the  highest  maximum  power  density 
of  152mWcm-2  (550 mAcm-2  @  0.277V),  whereas  the  FI -type 
flow  channel  (channel  width:  0.7  mm,  rib  width:  0.5  mm)  and  F3- 
type  flow  channel  (channel  width:  1.5  mm,  rib  width:  0.5  mm) 
had  a  maximum  current  density  of  580  and  540 mAcm-2  at  0.2  V 
and  a  maximum  power  density  of  132mWcm-2  (480 mAcm-2 
@  0.276V)  and  131mWcm-2  (450mA cm-2  @  0.289V),  respec¬ 
tively.  In  terms  of  the  constant  channel  width  (1 .0  mm),  the  F2-type 
flow  channel  with  a  narrow  rib  width  (channel  width:  1.0  mm, 
rib  width:  0.5  mm)  showed  much  better  performance  than  the 
broad  rib  width  case,  i.e.,  the  F4-type  flow  channel  (channel  width: 
1.0  mm,  rib  width:  1.0  mm).  The  experimental  trend  agrees  well 
with  the  previous  study  on  PEMFC  [12]  in  which  a  narrow  rib 
width  of  0.5  mm  brought  about  the  performance  improvement. 
These  results  show  the  strong  effect  of  the  flow  channel  geome¬ 
try  on  the  performance  of  DMFCs.  In  addition,  these  experimental 
trends  were  correctly  captured  by  the  present  DMFC  model.  To 
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Table  3 

Dh,  v  and  Re  values  for  calculation  of  the  theoretical  pressure  drop  (A  P)  of  the  different  flow  channel  geometries  calculated  by  Darcy-Weisbach  equation  (50)  according  to 
current  density  (stoichiometry  A,  =  3). 


Channel  name  DH  (mm)  v  and  Re  values  (v/Re)  according  to  current  density 


100  mA  cm-2 

200  mA  cm-2 

300  mA  cm-2 

400  mAcrrr2 

500  mAcrrr2 

FI 

0.82 

1.32/74 

2.60/146 

3.93/221 

5.24/295 

6.52/376 

F2 

1.0 

0.92/63 

1.82/124 

2.75/188 

3.67/251 

4.57/312 

F3 

1.2 

0.61/50 

1.21/99 

1.83/150 

2.44/200 

3.04/250 

F4 

1.0 

0.92/63 

1.82/124 

2.75/188 

3.67/251 

4.57/312 

Table  4 

Cell  properties  and  operating  conditions. 


Description 

Value 

Thickness  of  anode  GDL 

190  x  10~6  m 

Thickness  of  anode  CL 

30  x 10“6  m 

Thickness  of  cathode  GDL 

235  x  10-6m 

Thickness  of  cathode  CL 

30  x 10-6  m 

Thickness  of  membrane 

127  x 10-6  m 

Thickness  of  bipolar  plate 

2  x  10“3m 

Porosity  of  GDLs,  £Gdl 

0.7 

Porosity  of  CLs,  £Gl 

0.7 

Volume  fraction  of  ionomer  in  CLs,  £mc 

0.23 

Permeability  of  GDLs,  I<gdl 

1.0  x  10-12  m2 

Permeability  of  GDLs,  Ka 

1.0  x  10“12  m2 

Hydraulic  permeability  of  a  membrane,  Kmem 

5.0  x  10-19  m2 

Contact  angle  of  GDLs  and  CLs,  0 

92° 

Anode  stoichiometry 

2.5 

Cathode  stoichiometry 

3.0 

Cell  operating  temperature 

60  °C 

Anode/cathode  inlet  pressure 

Atmospheric 

Inlet  methanol  concentration 

1000  mol  m3 

Table  5 

Physical  and  transport  properties  in  the  individual  cell. 

Description 

Value 

Ref. 

Reference  exchange  current  of  anode,  air^a 

1.2  x  106  Am-3 

Reference  exchange  current  of  cathode,  air^c 

1.407  xlO3  Am-3 

Reference  oxygen  molar  concentration,  Co2,re/ 

7.68  mo/m-3 

Anodic  transfer  coefficients  for  MOR,  aa 

0.5 

Cathodic  transfer  coefficient  for  the  ORR,  ac 

0.8 

Activation  energy  for  MOR,  Eaa 

35570 J kg-1 

Wang  and  Wang  [28] 

Activation  energy  for  ORR,  Eac 

73200 J  kg-1 

Parthasarathy  et  al.  [29] 

Dry  membrane  density,  pmem 

2000  kg  ur3 

Equivalent  weight  of  electrolyte  in  the  membrane,  EW 

1.1  kg  mol-1 

Faraday  constant,  F 

96487  C  mol”1 

Universal  gas  constant,  Ru 

8.314J(molK)-1 

Surface  tension,  cr 

0.0625  Nm-1 

Liquid  water  viscosity,  p} 

3.5  x  10-4(Ns)m-2 

Diffusion  coefficient  of  methanol  in  liquid,  DlMeQH 

|Q-5.4163-999.778/Tm2  S~1 

Wang  and  Wang  [28] 

Diffusion  coefficient  of  methanol  in  gas,  D^e0H 

/ -6.954  x  10“2  + 4.5986  x  10“4T\  ,„_4  2  , 

(+9-4979  xlO-7!2  )  *  10  ™  S 

Yaws  [33] 

Diffusion  coefficient  of  oxygen  in  gas,  Dj, 

1.775  x  10-5(T/273)L823(1.0  1  3  x  105/p)m2  s-1 

Cussler  [34] 

Diffusion  coefficient  of  water  in  gas, 

2.56  x  10-5(T/307)2334(1.013  x  105/p)m2  s”1 

Cussler  [34] 

Electro-osmotic  drag  coefficient  of  water,  nd  w 

2.5 

Ren  et  al.  [35] 

elucidate  this  phenomenon,  it  should  be  noted  that  the  perfor¬ 
mance  of  the  cells  in  the  low  current  density  regions  was  hardly 
affected  by  the  channel  and  rib  width,  whereas  the  cell  performance 
in  the  high  current  density  regions  was  significantly  dominated  by 
the  channel/rib  cross  section  geometries.  This  could  be  ascribed  to 
the  differences  in  electrical  conduction  of  the  cells  with  different 
rib  areas,  which  influenced  the  ohmic  loss  in  the  IR  region  or  it 
could  be  caused  by  the  difference  in  the  mass  transport  capability 
of  the  flow  channel  geometries,  which  affected  the  access  of  the 
reactant  gas  (methanol  and  air)  to  the  electrode  layer  as  well  as 
the  product  removal  (H20  and  C02)  from  the  electrode  layer.  It  is, 
however,  very  difficult  to  elicit  a  conclusion  with  only  the  polar¬ 
ization  curves.  Thus,  the  numerical  and  experimental  studies  were 
chosen  to  accompany  the  polarization  curves  in  order  to  examine 


any  effects  of  the  flow  channel  geometries  on  the  performance  of 
DMFCs. 

Fig.  4  illustrates  the  simulated  current  density  distributions 
(Am-2)  over  the  membrane  at  the  operating  current  density  of 
400  mA  cm-2  with  F1-F4  type  channel  geometries.  These  results 
exhibit  a  wealth  of  current  non-uniformity  features  over  the  entire 
cell  of  9  cm2.  First,  it  is  found  that  the  current  density  distribu¬ 
tion  is  mainly  controlled  by  the  effect  of  oxygen  depletion  along 
the  cathode  downstream  rather  than  methanol  depletion  along  the 
anode  flow  path.  Second,  the  local  current  density  near  the  flow 
channels  is  higher  than  near  the  ribs,  implying  a  strong  rib  effect. 
Comparing  Fig.  4(b)  and  Fig.  4(d),  a  small  current  density  gradi¬ 
ent  between  the  channel  and  rib  regions  is  observed  in  the  F2-type 
flow  channel  (Fig.  4(b))  owing  to  its  narrower  rib  width  (0.5  mm). 
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Fig.  3.  Experimental  and  simulated  polarization  curves 
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(a)  F1-,  (b)  F2-,  (c)  F3-  and  (d)  F4-types  of  geometries. 
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Fig.  5.  Experimental  ohmic  resistance  of  the  cells  with  different  percentages  of  rib 
area  at  current  densities  of  100  mA cm-2  and  400  mA cm-2. 


On  the  other  hand,  the  broad  rib  of  the  flow  channel  (F4-type  flow 
channel)  created  a  little  longer  convection  route  and  thereby  led 
to  insufficient  sub-rib  convection  of  methanol  and  air.  This  result 
indicates  that  the  narrow  rib  of  the  flow  channel  is  advantageous 
to  the  two-phase  transport  under  the  ribs,  leading  to  the  increase 
of  the  cell  performance.  Based  on  the  same  rib  width  of  0.5  mm, 
a  comparison  of  Fig.  4(a),  Fig.  4(b)  and  Fig.  4(c)  indicates  that  the 
flow  channel  width  (Fl-type:  0.7mm,  F2-type:  1.0mm,  F3-type: 
1.5  mm)  is  another  factor  affecting  current  density  distribution  in 
the  membrane.  The  flow  channel  width  determines  the  direct  con¬ 
tact  area  between  the  fuels  and  the  GDLs,  and  the  broad  channel 
results  in  less  reactant  transport  loss  under  similar  active  electrode 
areas  (Fl-type:  885mm2,  F2-type:  840mm2,  F3-type:  870mm2). 
Therefore,  F3-type  with  the  broader  channel  width  shows  better 
current  density  distribution  in  the  flow  channel  than  Fl-type  and 
F2-type  with  the  narrower  channel  width.  Flowever,  due  to  the 
smaller  percentage  of  rib  area  in  F3-type,  i.e.  23.7%  and  the  resul¬ 
tant  higher  electronic  resistance,  the  F2-type  flow  channel  shows 
better  current  density  distribution  than  the  F3-type  flow  channel. 


Fig.  6.  Experimental  anode  and  cathode  polarization  curves  of  the  cells  using  the 
different  flow  channel  geometries  at  an  operating  temperature  of  60  °C. 

These  results  verify  that  the  narrow  rib  width  and  the  suitable  chan¬ 
nel/rib  aspect  ratio  (here,  the  channel/rib  aspect  ratio  is  2:1)  need 
to  increase  the  cell  performance  of  serpentine  DMFCs. 

4.2.  Parameters  of  the  flow  channel  geometry  on  performance  of 
DMFCs 

The  ohmic  resistance  of  these  cells  with  the  different  rib 
areas  was  experimentally  measured  at  a  low  current  density 
(100mA cm-2)  and  a  high  current  density  (400mA cm-2),  respec¬ 
tively,  as  presented  in  Fig.  5.  Increase  of  the  rib  area  resulted 
in  a  linear  decrease  of  the  ohmic  resistance  of  the  cells  (Fl- 
type:  360  mm2,  F2-type:  270mm2,  F3-type:  210mm2,  F4-type: 
420  mm2).  This  result  originates  from  the  difference  in  current 
passage  areas,  which  are  formed  by  contact  between  a  rib  and  a 
GDL.  It  should  be  noted,  however,  that  the  difference  in  the  ohmic 
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Fig.  8.  Experimental  and  simulated  pressure  drop  in  the  cathode  channel  of  F1-F4 
air  stoichiometry  X  =  3. 
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resistance  of  the  different  rib  dimensions  was  only  1.9  m£2  and 
that  the  voltage  loss  corresponding  to  the  ohmic  resistance  was 
merely  1. 7-6.8  mV  ( R0hmic =  VfossliA).  This  indicates  that  the  differ¬ 
ence  in  the  cell  resistance  generated  by  the  flow  channel  geometries 
less  seriously  affects  the  IR  loss  of  the  polarization  curve.  In  spite 
of  the  different  current  passage  areas  of  each  cell,  the  cell  per¬ 
formance  is  not  greatly  affected  by  the  ohmic  resistance  of  the 
cells  [37]. 

The  experimentally  measured  polarization  curves  of  the  anode 
and  cathode  shown  in  Fig.  6  describe  the  anode  and  cathode  reac¬ 
tion  features  of  the  cells  with  the  different  flow  channel  geometries. 
As  presented  in  Fig.  6,  the  performance  of  the  cells  with  these 
flow  channels  was  primarily  determined  by  the  cathode  overpo¬ 
tential  loss  rather  than  by  the  anode  overpotential  loss.  The  anode 
polarization  curves  of  all  the  cells  were  almost  the  same  in  the 
current  density  regions  below  500  mAcm-2,  whereas  the  cathode 
polarization  curves  differed  in  the  current  density  regions  above 
300  mAcm-2.  The  cathode  overpotential  loss  of  the  cells  with  F1-, 
F3-,  and  F4-type  flow  channels  gradually  increase  in  the  current 
density  regions  above  300 mAcm-2.  This  result  implies  that  the 
performance  difference  of  the  cells  with  these  flow  channel  geome¬ 
tries  could  primarily  be  ascribed  to  mass  transport  and  water 
removal  on  the  cathode,  since  convective  transport  governing  mass 
transport  depends  on  flow  channel  geometry  and  flow  conditions. 
Flowever,  it  is  observed  that  the  anode  overpotential  of  F3-type 
case  becomes  higher  than  others  above  about  500  mAcm-2,  which 
indicates  that  the  F3-type  flow  channel  is  not  likely  to  provide 
sufficient  methanol  demand  for  MOR  in  the  high  current  density 


Fig.  10.  Experimental  power  density  of  the  cells  according  to  the  stoichiometry  of 
supplied  air  in  the  cathode  at  an  operating  temperature  of  60  °C. 


region.  Although  an  increase  of  the  channel  width  would  expand 
the  direct  contact  areas  between  the  methanol  and  GDL,  a  broad 
channel  could  result  in  a  decrease  of  the  fuel  velocity  and  a  pres¬ 
sure  drop  within  the  channel,  and  it  finally  appears  to  limit  the 
convection  process  of  supplying  fuel  under  the  rib  and  to  remove 
the  produced  C02  from  the  cell. 
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Fig.  11.  Simulated  methanol  concentration  contours  (molirr3)  in  the  (a)  anode  flow  channel,  (b)  anode  GDL  and  (c)  anode  CL  at  current  density  of  400 mAcm-2  (F2-type 
channel). 
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Fig.  12.  Simulated  oxygen  concentration  contours  (mol  m-3)  in  the  (a)  cathode  flow  channel,  (b)  cathode  GDL  and  (c)  cathode  CL  at  current  density  of  400  mAcm-2  (F2-type 
channel). 


Nyquist  plots  of  the  impedance  spectra  for  anode  and  com¬ 
plete  cells  are  experimentally  measured  to  investigate  the 
characteristics  of  the  electrochemical  reaction  for  cells  using 
the  different  flow  channel  geometries.  Fig.  7(a)  and  (b)  shows 
the  Nyquist  plots  of  the  impedance  spectra  for  the  anode  at  a 
current  density  of  100 mAcm-2  and  400 mAcm-2,  respectively, 
and  Fig.  7(c)  and  (d)  for  complete  cells  at  a  current  density  of 
1 00  mA  cm-2  and  400  mA cm-2,  respectively.  The  impedance  spec¬ 
tra  of  the  anodes  (Fig.  7(a)  and  (b))  show  that  the  impedance  circles 
mostly  coincide  at  a  low  current  density  of  100 mAcm-2,  whereas 
the  Nyquist  plots  at  a  high  current  density  of  400  mAcm-2  are  def¬ 
initely  distinguished  at  low  frequencies  (below  1  Hz).  According  to 
Mueller  et  al.  [38],  the  low-frequency  arc  of  the  anode  impedance, 
which  can  be  observed  in  the  Nyquist  plots  at  frequencies  below 
1  Hz,  is  clearly  related  to  the  mass  transport  of  methanol  in  the 
anode.  Thus,  the  Nyquist  plots  at  a  high  current  density  confirm 
that  the  F3-type  flow  channel  has  a  larger  mass  transport  limitation 
than  the  other  flow  channel  geometries.  This  behavior  is  entirely 
consistent  with  the  result  of  the  anode  polarization  curve  of  Fig.  6. 
In  addition,  in  the  case  of  the  impedance  spectra  of  the  complete 
cells  (Fig.  7(c)  and  (d)),  the  mass  transport  limitation  from  FI -type 
and  F4-type  flow  channels,  as  well  as  from  the  F3-type  flow  chan¬ 
nel,  is  observed  at  a  high  current  density  of  400 mAcm-2.  These 
impedance  results  indicate  that  the  F2-type  flow  channel  exhibits 
better  mass  transport  capability  on  both  of  the  anode  and  cathode 
than  the  other  types. 

The  pressure  drop  of  the  cathode  channel  is  a  very  important 
factor  because  it  affects  the  two-phase  flow  associated  with  water 


removal  and  water  flooding  [39].  A  pressure  difference  between 
the  inlet  (pin)  and  the  outlet  (p0ut)  drives  the  fluid  flow.  Increasing 
the  pressure  drop  between  the  inlet  ( pin )  and  the  outlet  (p0ut)  will 
increase  the  mean  fuel  velocity  in  the  channel  and  improve  the 
convection.  However,  at  extremely  high  gas  velocities,  convective 
mixing  may  penetrate  into  the  electrode  itself,  causing  the  GDL  to 
retreat.  Thus,  convective  transport  that  dominates  mass  transport 
depends  on  the  flow  conditions  and  the  flow  channel  geometry 
[32].  According  to  Barbir  et  al.  [40],  pressure  drop  can  be  used  as  a 
diagnostic  tool  to  detect  flooding  in  a  fuel  cell.  The  pressure  drop 
in  the  cathode  channel  increases  with  the  current  density.  This  is 
mainly  due  to  either  liquid  water  droplets  or  water  films  in  the  flow 
channels,  which  are  related  to  the  amount  of  water  that  is  produced 
from  the  oxygen  reduction  reaction  based  on  Faraday’s  law  [41 ,42]. 
In  addition  to  this  general  fact,  it  should  be  noted  that  the  cathode 
pressure  drop  could  be  affected  by  the  flow  channel  geometry. 

Fig.  8  shows  the  experimentally  measured,  analytically  calcu¬ 
lated,  and  numerically  simulated  pressure  drops  along  the  cathode 
flow  path  of  F1-F4  type  geometries  at  the  operating  current  densi¬ 
ties  from  1 00  to  500  mA  cm-2  under  air  stoichiometry  of  A  =  3.  The 
analytical  pressure  drop  in  the  flow  channels  of  the  cathode  was 
calculated  by  the  Darcy-Weisbach  equation  (51),  which  enables 
the  calculation  of  the  pressure  loss  due  to  friction  within  a  given 
length  of  the  channel.  The  experimental  pressure  drop  was  directly 
measured  from  the  single  cells  with  the  different  flow  channels. 
From  Darcy-Weisbach  equation  (51),  the  analytical  pressure  drop 
in  the  flow  channels  is  linearly  proportional  to  the  velocity  of  the 
flow  and  the  channel  length.  This  is  illustrated  well  in  Table  3  and 
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Fig.  13.  Simulated  liquid  saturation  contours  in  the  (a)  anode  GDL,  (b)  anode  CL  (c)  cathode  CL  and  (d)  cathode  GDL  (F2-type  channel). 


Fig.  8.  Furthermore,  the  simulated  results  of  Figs.  8  and  9,  which 
show  the  contour  of  pressure  drop  in  the  cathode  channel  of  F1-F4 
type  geometries  at  an  operating  current  density  of  400  mA  cm-2, 
exhibit  an  analogous  tendency  on  the  cathode  pressure.  Flowever, 
the  calculated  and  simulated  pressure  drops  are  smaller  than  the 
pressure  drop  experimentally  measured  using  a  nanometer,  which 
indicates  that  the  water  flooding  in  the  cathode  channel  is  signifi¬ 
cant  in  the  experimental  cells.  As  shown  in  Fig.  8,  the  pressure  drop 
in  the  F2-type  flow  channel  is  somewhat  lower  than  the  pressure 
drop  of  the  F4-type  flow  channel,  which  has  the  same  flow  velocity 
but  a  shorter  channel  length  than  the  F2-type  flow  channel.  Gener¬ 
ally,  during  normal  operation,  liquid  water  initially  accumulates  in 
the  GDL  under  the  ribs,  as  ribs  are  the  coolest  location  [43].  Yamada 
et  al.  [44]  found  that  liquid  water  accumulated  more  rapidly  under 
ribs  at  a  higher  current  density  and  the  liquid  water  under  ribs  was 
more  difficult  to  expel  due  to  its  longer  transport  route.  Therefore, 
considering  the  relationship  between  the  pressure  drop  and  the 
two-phase  flow  in  the  cathode,  this  result  implies  that  an  F2-type 
flow  channel  with  a  narrow  rib  eliminates  the  accumulated  liq¬ 
uid  water  under  the  ribs  more  rapidly.  In  addition,  the  good  water 
removal  capability  provides  shorter  convection  routes  to  supply 
air  quickly  at  the  sites  from  which  the  water  was  removed,  which 


improves  the  mass  transport  characteristics  in  the  high  current- 
density  region.  These  characteristics  of  the  F2-type  flow  channel 
enable  the  F2-type  channel  to  have  higher  cell  performance  than 
the  other  geometric  channels.  On  the  other  hand,  the  higher  pres¬ 
sure  drop  of  the  FI -type  flow  channel  could  be  caused  by  the  long 
water  column  from  the  narrow  channel;  the  higher  pressure  drop 
of  the  F4-type  flow  channel  could  result  from  the  liquid  water  that 
accumulated  under  the  ribs  and  from  the  longer  convection  routes 
due  to  the  broad  rib.  These  phenomena  appear  to  make  it  much 
more  difficult  for  the  two  channels  to  expel  the  accumulated  liq¬ 
uid  water  outside  the  cells.  The  F3-type  flow  channel  experienced 
a  lower  pressure  drop  in  the  channel  due  to  a  broad  channel  width 
and  a  slow  flow  velocity,  making  it  more  likely  to  have  a  dead  site 
under  the  ribs.  Thus,  the  F2-type  flow  channel  geometry  is  bene¬ 
ficial  for  the  operation  of  the  cell  since  because  it  easily  removes 
excess  liquid  water  from  the  cell  and  contributes  to  mitigating  the 
mass  transport  problem  in  a  higher  current  density. 

Effects  of  the  stoichiometry  of  air  supplied  in  the  cathode  on 
the  experimentally  measured  power  densities  of  the  cells  are 
presented  in  Fig.  10.  Generally,  excess  air  should  be  supplied  to 
the  cathode  to  sustain  the  electrochemical  reaction  and  to  prevent 
oxygen  starvation;  however,  this  can  influence  the  two-phase  flow, 
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which  depends  on  the  pressure  drop,  fuel  distribution,  and  water 
flooding.  Therefore,  the  effects  of  the  stoichiometry  of  air  in  the 
cathode  can  vary  with  the  flow  channel  geometries.  As  shown  in 
Fig.  10,  the  performance  of  the  cells  with  F1-,  F2-,  and  F4-type  flow 
channels  was  improved  with  an  increase  of  the  stoichiometry  of 
the  supplied  air  in  the  cathode.  These  results  show  that  increasing 
the  stoichiometry  of  the  air  in  the  cathode  leads  to  better  fuel  dis¬ 
tribution  and  more  efficient  removal  capability  of  product  water.  In 
contrast,  the  performance  of  the  cell  with  an  F3-type  flow  channel 
was  not  considerably  enhanced  by  increasing  the  air  stoichiome¬ 
try.  It  is  believed  that  the  lower  pressure  drop  in  the  flow  channel, 
caused  by  the  broad  channel  width,  leads  to  difficulty  in  the  supply 
of  the  fuel  under  the  rib  and  in  the  removal  of  the  produced  water 
from  the  cell  despite  the  increase  in  the  amount  of  the  air. 

4.3.  Detailed  methanol,  air,  liquid  saturation  contours  for  the 
F2-type  flow  channel 

Besides  the  comparison  of  the  predicted  results  and  measured 
data,  multidimensional  distributions  of  methanol,  oxygen,  and  liq¬ 
uid  saturation  are  also  obtainable  from  the  simulations  and  should 
provide  greater  insight  into  the  operation  of  DMFCs.  Fig.  1 1  shows 
the  numerically  simulated  methanol  concentration  (molm-3)  on 
the  anode  flow  channel  (Fig.  11(a)),  anode  GDL  (Fig.  11(b))  and 
anode  CL  (Fig.  11(c)).  The  methanol  concentration  was  calcu¬ 
lated  with  F2-type  channel  geometry  (channel  width:  1.0  mm,  rib 
width:  0.5  mm)  at  a  current  density  of  400  mA cm-2.  In  compar¬ 
ison,  the  methanol  concentration  in  DMFCs  decreases  from  the 
anode  flow  channel  (Fig.  11(a))  toward  anode  CL  (Fig.  11(c))  along 
the  through-plane  direction.  In  particular,  the  steep  gradient  of 
methanol  concentration  between  the  anode  GDL  and  the  anode  CL 
is  observed  and  is  due  to  a  limited  diffusion  rate  of  methanol  from 
the  anode  flow  channel  into  the  anode  CL  as  well  as  fast  consump¬ 
tion  of  methanol  at  the  CL.  The  limited  diffusion  rate  of  methanol 
can  be  ascribed  to  presence  of  the  C02  gas  in  the  flow  channel  and 
GDL.  According  to  Lu  et  al.  [45],  the  C02  gas  bubbles  are  held  on  the 
carbon  paper  by  strong  surface  tension  until  they  grow  into  larger 
slugs  for  detachment.  Once  the  bubbles  grow  to  a  sufficient  size, 
they  detach  and  sweep  along  the  GDL  surface  of  the  channel.  FIow- 
ever,  when  C02  gas  cannot  be  removed  efficiently  from  the  surface 
of  the  backing  layer  and  they  remain,  the  GDL  and  flow  channel 
blocked  by  C02  gas  bubbles  and  slugs  can  restrict  the  through-plane 
transport  of  methanol  to  the  anode  CL,  thereby  increasing  the  con¬ 
centration  polarization  [46].  This  bubble  motion  is  governed  by  the 
momentum  of  liquid  flow  as  well  as  the  surface  adhesion  between 
bubbles  and  substrate  [45].  Thus,  the  anode  flow  channel  can  be  an 
important  factor  to  obtain  better  mass  transport  of  methanol  and 
C02  gas  on  the  DMFC  anode.  The  importance  of  the  flow  channel 
was  also  confirmed  from  the  in-plane  transport  of  methanol.  From 
Fig.  11(b)  and  (c),  the  methanol  concentration  decreased  near  the 
methanol  outlet  and  around  the  channel  ribs.  These  phenomena  are 
affected  by  the  sub-rib  convection  of  the  serpentine  flow  channel, 
the  flow  velocity  and  the  pressure  drop.  This  indicates  that  the  mass 
transport  of  methanol  and  C02  gas  on  the  DMFC  anode  depends  on 
the  dynamic  mass  transfer  of  the  anode  flow  channel. 

Fig.  12  exhibits  the  numerically  simulated  air  concentration 
(molm-3)  on  the  cathode  flow  channel  (Fig.  12(a)),  cathode  GDL 
(Fig.  1 2(b))  and  cathode  CL  (Fig.  1 2(c))  with  F2-type  channel  geom¬ 
etry.  The  cell  condition  is  the  same  as  that  of  Fig.  11.  It  is  seen  in 
Fig.  12  that  the  oxygen  concentration  gradient  along  the  through- 
plane  direction  is  smaller  than  the  methanol  concentration  (Fig.  11), 
mainly  because  of  the  larger  gas-phase  diffusivity  of  oxygen.  Fig.  13 
displays  the  liquid  saturation  contours  in  both  anode  and  cathode 
GDLs  and  CLs  for  the  F2-type  channel  geometry.  It  is  clearly  seen 
that  the  liquid  phase  is  dominant  in  the  anode  porous  region  (GDL 
and  CL)  whereas  the  cathode  GDL  and  CL  are  gas-dominated.  The 
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Fig.  14.  Experimental  long-term  stability  (Fig.  13(a))  and  polarization  curves 
(Fig.  13(b))  of  a  large  cell  (150  cm2)  with  the  F2-type  flow  channel  geometry. 

fact  further  indicates  that  C02  in  the  anode  and  water  in  the  cathode 
are  efficiently  removed  with  the  F2-type  channel  geometry. 


4.4.  Long-term  stability  test  ofF2-type  flow  channel 

The  experimental  results  of  long-term  stability  and  the  polariza¬ 
tion  curve  of  a  large  cell  (150  cm2)  with  the  F2-type  flow  channel 
are  given  in  Fig.  14(a)  and  (b),  respectively.  The  long-term  test 
of  the  cell  was  performed  with  a  constant  current  mode  of  20  A 
(133  mA cm-2)  at  a  methanol  flow  rate  of  A,  =  2.5  and  an  air  flow 
rate  of  A  =  3.0  for  700  h.  As  shown  in  Fig.  14(a),  the  large  cell 
with  the  F2-type  flow  channel  exhibited  a  constant  performance 
of  9.60 W  (0.480V,  64.0 mW cm-2)  for  700 h.  At  the  beginning 
of  the  test,  the  cell  provided  output  power  of  9.66  W  (0.483  V, 
64.4 mW cm-2),  and  this  output  was  maintained  for  700 h.  The 
observed  performance  loss  over  the  700  h  period  was  scarcely 
(9.66  W  (0.483  V)  9.60  W  (0.48  V)).  This  long-term  stability  of 
the  cell  was  also  confirmed  from  the  polarization  curves  of  the 
cell  before  and  after  the  long-term  test.  As  given  in  Fig.  14(b), 
the  polarization  curves  of  the  cell  before  and  after  the  long¬ 
term  test  for  700  h  were  almost  the  same:  the  performance  loss 
for  a  nominal  voltage  of  0.40  V  was  scarcely  (15.65W  (39.13  A, 
1 04.33  mW  cm-2 )  -*  1 5.62  W  (39.04  A,  1 04.1 3  mW  cm”2))  after  the 
long-term  test  of  700  h.  These  results  indicate  that  the  cell  with  the 
F2-type  flow  channel  provides  good  long-term  stability  as  well  as 
high  performance  (0.4  V,  104.13  mW  cm-2). 

According  to  the  previous  paper  [31,47],  performance  loss  of 
fuel  cells  in  long-term  operation  mainly  originates  from  cathode 
degradation  caused  by  electrocatalyst  agglomeration  or  deteriora¬ 
tion  of  mass  transport  within  the  GDLs  and  the  CLs.  Based  on  these 
suggestions,  our  results  suggest  that  the  F2-type  flow  channel  pro¬ 
vides  very  good  fuel  distribution,  mass  transport  and  removal  of 
the  products  (H20  and  C02)  in  DMFCs  during  long-term  operation. 


5.  Conclusion 


In  this  study,  the  effect  of  flow  channel  geometry  on  the  perfor¬ 
mance  of  DMFCs  was  numerically  and  experimentally  investigated 
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using  single-path  serpentine  patterns  with  different  channel 
widths,  rib  widths,  and  channel  lengths.  In  spite  of  the  different 
current  passage  areas  of  each  cell,  the  cell  performance  was  not 
affected  by  the  resistance  of  the  cells.  The  anode/cathode  polar¬ 
ization  curves  showed  that  the  performance  of  the  cells  with  the 
different  flow  channel  geometries  was  primarily  determined  by  the 
cathode  overpotential  loss  rather  than  by  the  anode  overpotential 
loss,  particularly  in  the  current  density  regions  above  300  mA  cm-2. 
According  to  impedance  spectra,  the  F2-type  flow  channel  had  bet¬ 
ter  mass  transport  capability  in  both  the  anode  and  cathode  sides 
than  the  other  types.  The  pressure  drops  measured  in  the  cath¬ 
ode  channel  indicated  that  the  F2-type  flow  channel  eliminated 
liquid  water  accumulated  under  the  ribs  more  readily;  the  good 
water  removal  capability  also  provided  short  convection  routes 
to  supply  air  quickly  at  the  removed  water  sites.  Finally,  such 
advantages  resulted  in  mass  transport  improvement  of  the  F2-type 
flow  channel  in  high  current  density  regions.  In  addition,  the  long¬ 
term  operation  test  of  a  large  cell  (150  cm2)  with  the  F2-type  flow 
channel  showed  that  performance  degradation  of  the  cell  hardly 
occurred  for  700  h.  This  result  confirmed  that  F2-type  flow  channel 
provided  very  good  fuel  distribution,  mass  transport  and  removal  of 
the  products  (H20  and  C02)  in  a  DMFC  during  long-term  operation. 
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